Asymmetry between the electron- and hole-doped Mott 
transition in the periodic Anderson model 
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We study the doping driven Mott metal-insulator transition (MIT) in the periodic Anderson model 
set in the Mott-Hubbard regime. A striking asymmetry for electron or hole driven transitions is 
found. The electron doped MIT at larger U is similar to the one found in the single band Hubbard 
model, with a first order character due to coexistence of solutions. The hole doped MIT, in contrast, 
is second order and can be described as the derealization of Zhang-Rice singlets. 

PACS numbers: 71.30.+h, 71. lO.Fd, 71.27. +a 



I. INTRODUCTION 

The nature of the Mott transition, i.e. the metal- 
insulator transition driven by electronic correlations, is 
a central problem in physics of strongly correlated elec- 
trons systems. The relevance of the problem was initially 
emphasized by Mott 1 in the 40's, trying to explain why 
some materials with odd electrons per the unit cell, like 
NiO, are insulators. In a Mott metal-insulator transi- 
tion (MIT), a metallic system with a partially filled elec- 
tron band suddenly opens an insulating gap. In practice 
the transition is usually driven by temperature, applied 
pressure, or chemical doping. The origin of the mecha- 
nism is in the correlation effects due to the strong on-site 
Coulomb repulsion experienced by electrons occupying 
rather localized orbitals, such as d in transition metal 
oxides or / in heavy fermion compounds. 

The classical example of an experimental system ex- 
hibiting a Mott transition is vanadium oxide V2O3, which 
has received continuous attention since the pioneering 
work of McWhan 2 in the 70's. That compound has a 
finite temperature first order metal-insulator transition 
that terminates at a high temperature 2nd order critical 
point, in analogy with the finite liquid-gas transition line 
in water. 

From a theoretical perspective, it is considered that 
the Hubbard model, which contains a tight binding band 
plus an interaction term that describes local Coulomb re- 
pulsion between electrons occupying the same site, is a 
minimal Hamiltonian that may capture the basic physics 
of the Mott MIT. The most significant work on this model 
was initiated by Hubbard 3 in the 60's, where, starting 
from the insulating state at large interaction values, he 
described how the system may close the correlation gap 
as the bandwidth is increased to values of the order of the 
Coulomb repulsion energy. Later, in the 70's Brinkman 
and Rice , using a variational approach, started from the 
metallic state and showed how it can be destroyed by in- 
creasing correlation effects when the interaction strength 
becomes of the order of the bandwidth. Finally in the 
90 's the theoretical development of the dynamical mean- 



field theory (DMFT) 5 ' 6 allowed to get new insight on this 
problem. In the scenario for the Mott transition realized 
in the DMFT solution of the Hubbard model, for low 
temperatures and moderate interaction, the half filled 
Mott insulator can be driven to a correlated metallic state 
through a first order transition 6 . This transition can oc- 
cur as a function of correlation strength, temperature 
or doping. The solution of the Hubbard model within 
DMFT provided not only a connection between the ap- 
proaches of Hubbard and Brinkman- Rice by showing how 
the system evolves from a metal to an insulator, but also 
produced a detailed description of the basic experimental 
phenomenology observed in the V2O3 compound 6 . In ad- 
dition, it was later shown that the MIT can be described 
in terms of a Ginzburg-Landau scenario 7,8,9 with theo- 
retical predictions for the critical behavior of obscrvables 
near the 2nd order critical point, that were eventually 
confirmed by experiments 10 . 

The physics of the Mott transitions, especially 
those driven by doping at low temperatures, became 
of unparalleled interest in condensed matter physics 
with the discovery of the high temperature cupratc 
superconductors 11 in the 80's and, in smaller but also sig- 
nificant measure, by the discovery of the non-Fermi liquid 
behavior in heavy fermion systems 12 . In those systems, 
the effect of strong correlations is undisputed, since the 
active electronic degrees of freedom involve the localized 
d and / orbitals. Therefore, these systems are identified 
as doped Mott insulators, however their phase diagrams 
and the evolution of their physical properties cannot be 
associated to the DMFT scenario for the Mott transition 
that was so successfully applied to V2O3. 

Besides the Hubbard model, the periodic Anderson 
model (PAM) is another minimal Hamiltonian that is 
often investigated in the context of strongly correlated 
electron systems. That model contains two orbitals per 
unit cell, one local with on-site Coulomb repulsion and 
the other non-interacting and itinerant. At each lattice 
site, the two orbitals are hybridized. This model is more 
realistic than the Hubbard, since it describes with greater 
detail the actual situation in real compounds. For in- 



stance, in transition metal oxides where the overlaps be- 
tween neighboring oxygen p— orbitals provide itineracy to 
the electrons, while the localized d— orbitals of the tran- 
sition metal experience the stronger correlation effects. 

Despite the higher degree of detail included in the 
PAM, it is often assumed in theoretical approaches that 
the physical behavior of the PAM would result qualita- 
tively similar to that of the Hubbard model at low fre- 
quencies. That statement can be mathematically jus- 
tified in certain parameter regimes, however its general 
validity is less evident 

The DMFT is a theoretical approach that is mathemat- 
ically exact in the limit of large lattice dimensionality 5 ' 6 
which has been extensively used to study the Mott tran- 
sition in the Hubbard model, and, to a lesser degree, 
has also been employed to investigate the physics of the 
PAM 15,16,17 ' 18 . Therefore, in the light of the previous dis- 
cussion, a natural question to address is whether within 
DMFT the Mott transition scenario of the PAM is indeed 
qualitatively similar to that of the Hubbard model or, if 
contrary to usual expectations, it brings about new phys- 
ical behaviors. This issue has been the focus of our recent 
investigations 19 , where we showed that in fact a different 
type of doping driven Mott MIT is realized in the PAM, 
even in a parameter regime whether it might be expected 
that the identification with the Hubbard model may hold. 
The present study extends and provides further details to 
that work. In particular, we present new data that illus- 
trate the different behavior of physical quantities on the 
two sides of the transition, we provide new comparisons 
of the Green functions to results using a T = numerical 
technique, and we extend the discussion of the physical 
origin of two transitions. We should clarify that there is a 
qualitatively different insulator state that can be realized 
in the PAM, namely, the Kondo insulator. It is obtained 
in the particle-hole symmetric case, where the total oc- 
cupation is even (ie, 2 electrons per site). The Kondo 
insulator is qualitatively different from a Mott insulator, 
because it is due to band hybridization effect and cor- 
relations merely serve to reduce the hybridization gap. 
This insulator has a temperature-driven crossover to a 
metallic state that has been investigated within DMFT 
using QMC techniques 13,14 . In contrast, the Mott insu- 
lator state that we consider here corresponds to a state 
with odd total number of particles and its physical origin 
is entirely due to strong correlation effects. 

The paper is organized as follows. In Sec. II we in- 
troduce the PAM and justify the choice of the param- 
eter regime. We also summarize the DMFT equations 
and provide details on the numerical techniques we use 
to solve the associated impurity problem. In Sec. Ill 
we present the results and discuss the Mott transitions 
found in the PAM. In Sec. IV we present a discussion of 
the physical origin of the different scenario for the MIT 
found in the PAM with respect to the HM. In Sec.V we 
present the conclusions. 
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FIG. 1: Schematic representation of the periodic Anderson 
model for the case of a one dimensional chain. At each site of 
the lattice, a d orbital (open circle) hybridizes with a p orbital 
(full circle) through the amplitude tpd (dot-dashed line) . Two 
electrons in the d orbital experience a Coulomb repulsion U. 
The hopping amplitude between the p orbitals at neighboring 
sites i and j is ty = t (solid line). 

II. METHODOLOGY 
A. Model 

The Hamiltonian of the PAM is given by 

H = ~ ^(ptxPjv + PjvPiv) + ( e P - ti^PiaPio 

<ij>cr icr 

+ ( £ d - M) d ^ dllJ + ^ ( d tP™ + Pt d i") 
icr icr 

i 

(i) 

Here pi„ and pf a operators destroy and create electrons 
at p orbitals on site i with spin a. The p orbitals have 
site energy e p and overlap via the hopping term iy = t 
to form a band. di a and d\ a operators destroy and create 
electrons at d orbitals on site i with spin a. The local 
d orbitals have site energy and are hybridized to the 
p orbitals with the constant on-site amplitude tpd- U is 
the energy cost of double occupation of the d orbitals at 
each site and [i is the chemical potential. Ao = — e p , 
the difference between the d level and the center of the 
p band, defines the bare charge-transfer energy. Fig. 1 
shows a schematic representation of the Hamiltonian. 

The insulator solutions of the PAM have been stud- 
ied in literature mostly in the symmetric regime A = 0, 
fi = 0, where the system is a Kondo (or renormalized- 
band) insulator. In the present work, we shall focus on a 
different and less explored parameter regime, where the 
system is in a Mott insulator state. A key difference be- 
tween these insulator states is that the former is realized 
at an even total electron occupation (n tot = 2), while the 
latter is realized at an odd occupation (n tot = 1 or n to t = 
3). Yet, the Mott insulator state may still be classified 
as a charge-transfer insulator or Mott-Hubbard insulator, 
according to the Zaancn-Sawatzky- Allen scheme 20 . 
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The action associated with the Hamiltonian (1) reads 
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where ip ko - = {d a ,Pko-}, ip^o- = { d tiP kcr } and tnc inverse 
matrix propagator Gq 1 is given by: 

G^{k,iu Jn )=( lUJn+ J 1 - €d . *** V 
0<T V ' ' \ t pd icun + V - e p - e k J ' 

(3) 

where e k is the Fourier transform of the hopping term t. 
The lattice Green's function G is then written using the 
Dyson equation G _1 = Gq 1 — S, where 



T, a (k,iuj n ) 




(4) 



and Ti a {k 1 ioj n ) is the d— electron self-energy. Here we 
arc interested in a magnetically disordered state, thus 
the spin index can be dropped. In this case, the lattice 
Green's functions for the p and d electrons are explicitly 
given by: 



G \p}(k, ioj n ) =iu> n + fi - e p - e k 



pd 



(5) 



iuj n + /i - e d - £(fc, iuj n ) 



G d l(k, iuj n ) =-iio n + jj, - e d - E(fc, iuj n ) 



pd 



(G) 



iu n + (i - t p - e k 



The local Green's functions arc then obtained perform- 
ing the integration over momenta, 

G a (iu) n ) = ^-^2G a (k,iuj n ) = j p (e)G a (e,iuj n )de, 

k ^ 

(7) 

where a = pp,dd and Po(e) = ^2 k S(e — e k ) is the free 
(U = and t pd = 0) density of states of the p electrons. 

The PAM has some simply solvable limits such as of 
vanishing hybridization, t pd = 0, or of vanishing corre- 
lation strength, U = 0. For the latter case, the PAM 
describes two hybridized one-particle bands, that are ob- 
tained diagonalizing the Hamiltonian and read: 



E±(k) = | ( e d + e p + e k - 2^ ± x / {e k - A ) 2 + 4tJ 



(8) 



From the many interesting parameter regimes that this 
model has, we shall focus our study on a particular one 
where the low energy physics of the PAM is a priori 
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FIG. 2: Top left panel: solid lines are the two branches E±(e) 
for U — 0, Ao = 1, tpd — 0.9 and jj, = 0.529, that gives n to t = 
3. Top right panel: density of states for the p and d electrons 
(dashed and solid line), for the same model parameters as 
in the left panel. Note that the effective distance between 
the bands, A ~ 2, results to be larger than the bare value 
A = e d - e p = 1. The free ([7 = and t pd = 0) DOS of 
the conduction electrons is semi-elliptical with with a half- 
bandwidth equal to unity. Bottom panel: particle occupation 
rid (solid line), n p (dashed) and ntot (dotted) as a function of 
the chemical potential for [7 = 0, Aq = 1, t p d = 0.9. 



expected to correspond to that of the Hubbard model. 
Thus, we consider the case where the localized d orbital 
is near the Fermi energy and with an occupation close to 
one, while the p orbital band is well beneath in energy 
and almost fully occupied. By virtue of the hybridiza- 
tion term, the d electrons acquire a finite dispersion and 
form a narrow band that crosses the Fermi energy. It has 
a bandwidth ~ t^ d /A, where A is the distance between 
the two hybridized bands, A w E + — E- > Aq, and 
is subject to strong correlation effects when the on-site 
Coulomb term is turned on. 

For reference, the solution for the non-interacting case 
in the chosen parameter regime is shown in Fig. 2. We set 
the units adopting a model semicircular density of states 
for the p-electrons with half-bandwidth D = 2t = 1 As 
we shall see in the next section, this density of states is 
actually realized in a Bethe lattice in the limit of infinite 
spatial dimensions. In the upper panels we show the 
single particle dispersion E± and the resulting density 
of states for Ao = 1, t pd = 0.9 and /i = 0.529, that 
gives a total occupation per site ntot equal to 3. In the 
lower panel we plot the total occupation n to t = np + n d 
together with the partial occupation. The plateaux in 
the curves signal incompressible states that correspond to 
insulators. These are observed at n to t = 2 and n tot = 4. 
In the first case, ntot equal to 2, one has a hybridization 
gap insulator, since the band £L is full and E + is empty. 
For ntot — 4 the state is of a full band insulator. For 
2 < n to t < 4 the system is metallic, since the lower band 
with mostly p character is full and the narrow band with 
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mainly d character is partially filled. 

We shall now focus on the effects of correlation on such 
a metallic case, realized for n to t = 3, where the lower 
electron band with mostly p character is close to being 
full (i.e. is occupied by two electrons) and the narrow 
d electron band is half-filled (i.e. with an occupation rid 
close to 1). In the chosen parameter regime, there is a 
single band crossing the Fermi energy which has mostly 
d electron character. So correlations effects will affect it 
strongerly. For values of the interaction U larger than its 
bandwidth ~ tz^j A one may expect that a correlation 
gap would open and the system becomes a Mott insula- 
tor. This is the regime where the identification of the low 
energy physics of the PAM with the one band Hubbard 
model may hold. 



B. DMFT and the limit of infinite dimensions 

To go beyond this qualitative discussion we need to 
obtain reliable solutions of the model Hamiltonian in the 
strongly interacting regime. We thus recur to the DMFT 
formulation where exact numerical methods can be used 
to solve the problem 6 . 

The DMFT solution becomes exact in the limit of large 
spatial dimensionality or, cquivalcntly, large lattice con- 
nectivity z. For this limit to remain physical one is re- 
quired to rescale the hopping t amplitude as t/^fz, so 
that the density of states po{e) = J2k — £fc ) gi yes a S~ 
nite value for the mean kinetic energy 5 . As is well known 
in DMFT, the specific lattice structure is not essential, 
and several lattice types could be used. For instance, the 
free (i.e., t p d = 0, U = 0) density of states of the hyper- 
cubic lattice 5 ' 15 , which is the generalization of the square 
lattice to the limit of high z, reads: 

#■"<«> -iTf "*(-£)■ w 

where e denotes the nonintcracting single particle energy. 
Another lattice type which is often adopted is the Bcthe 
lattice, whose density of states reads, 



Po(e) = ^V^^- (10) 

In the following we shall adopt this type of lattice struc- 
ture, as it is better suited for some of the numerical meth- 
ods that we shall employ. As unit of energy we set the 
half-bandwidth of the Bcthe lattice semi-circular DOS, 
D = 2t = 1 . The key mathematical simplification arising 
from the z — ► oo limit is the locality of the self-energy, 
i.e. its k independence. Thus, there is no longer need to 
keep the momentum label in the single particle energies 
of the band structure e^, and the energy e itself is simply 
kept as the quantum number. 



C. Mean-field equations 

In the limit of large lattice connectivity z — * oo, the 
PAM can be exactly mapped onto a single impurity An- 
derson model supplemented with a self consistency condi- 
tion. The derivation of the DMFT equations has already 
been presented in detail elsewhere 6,21,22 , so here we shall 
just briefly summarize the main steps and the final ex- 
pressions. 

A direct way to derive the DMFT equations is to ap- 
ply the cavity method 6 . The key idea is to focus on a 
given (any) site of the lattice and to integrate out the 
degrees of freedom on all the other lattice sites in order 
to obtain the local effective action at the selected site. In 
doing that, one shall also obtain a self-consistency condi- 
tion which restores the translational invariance that was 
(temporarily) broken with the selection of a given lat- 
tice site. After integration of all sites other than the one 
selected, the local effective action is obtained, 

S off =- f dr [ ^ ^r'y;V+(T)a - 1 (r-r , )Va(r , ) 
Jo Jo a 

+ U / dr [n dT (r) - \] [n di (r) - ±] , 
Jo 

(11) 

where ip a = {dg-.p^}, ip+ = {<i+,p+} correspond to the 
two atomic orbitals of the given (arbitrary) site of the 
lattice. The local inverse propagator reads 

g-if iw \ =( iuj n + Li-£d tpd _ \ 

° \ tpd + /! — € p — t 2 G pp J 

(12) 

where G pp is the cavity Green's function that encodes 
the information of the propagation of electrons in the 
lattice, restricted not to return to the local site. These 
two equations define the so called associated impurity 
problem of the model. 

In order to restore the translational invariance and to 
obtain a closed set of equations, one has to relate the 
local inverse propagator Qo to the Green's function of 
the original lattice. In the Bethe lattice this relation 
is simple 6 : the p— electron cavity Green's function G pp 
becomes the local p— electron Green's function, G pp = 
G pp . From this relation and (12), one obtains the self- 
consistency condition for the impurity problem. It can 
be casted only in terms of [^o]drf and reads, 

[Go^ddiiun) =iu n + M - e d 



iui n + [i — 6p — t 2 Gpp 

In practice an iterative procedure is implemented to 
solve the set of DMFT equations: given an ansatz for 
[SoW, and the fact that the interactions are local and 
only act on the d orbital, the impurity many-body prob- 
lem (11) can be solved to produce a local d— electron 
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Green's function G 
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This defines a lo- 



G d< 7, that allows for the 



cal self-energy £ = [Q Q 1 
calculation of the local p— electron Green's function G pp 
(6) 



Gp P (ioj n ) 



Po(e)de 



+ /- 1 e P iu) n + l i-e d -T,{iu> n ) 



(14) 

The obtained G pp and S are then used as input to the 
self-consistency condition Eq. 13 to produce a new [<?o]dd- 
This process is iterated until convergence is reached. At 
the self- consistent point, the Green's functions Gd d , G pp 
and the self-energy £ correspond to the local propaga- 
tors of the original lattice model. Moreover, since in the 
z — > oo limit S is local, then all the k dependent (or e 
dependent) propagators of the original lattice can also be 
computed from this self-energy. 

It is useful rewrite the local Green's functions in terms 
of the Hilbert transform of the density of states D(£) = 
depo(e) / — e) . For the conduction electron Green's 
function G pp we get 



Gp P (iuJn) =D 



iui n 



pd 



(15) 

and for the d— electron Green's function Gdd we obtain 

1 



Gdd(iu n ) 



iuj n + p - e d - E(iw„) 

tpd 

ito n + /i - e d - Tj{iuj n ) 



Gpp(iu! n ). 



(16) 



This expression has a transparent physical interpretation: 
there are two processes that a d— electron can undergo: 
either remain fluctuating at the local site (first term), 
or fluctuate for some time, then jump to the p site and 
propagate, and then return and fluctuate some more time 
(second term). 



D. Numerical solution 

For the solution of the associated impurity many-body 
problem (11), one may use a variety of techniques 6 . Here 
we shall employ two numerical methods which are both a 
priori exact: Hirsch-Fye quantum Monte Carlo (QMC) 23 
and exact diagonalization (ED). The interest of using dif- 
ferent techniques is that they have complementary range 
of applicability and that they allow for a crosscheck of 
the numerical results. The first method, QMC, is a finite 
temperature calculation and is exact in the statistical 
sense. The other method is formulated at T = and rely 
in a finite size representation of the local site environ- 
ment (i.e., the cavity Green's function) by a bath of non- 
interacting atomic sites connected to the local impurity. 



In the limit of large number of atomic sites in the bath, 
this approach also becomes a priori exact. Both methods 
have already been well documented in the literature 6 , so 
here we shall limit ourselves to briefly provide the rele- 
vant technical details. 

To implement the QMC, it is useful to first perform 
the integration on the non-interacting local p— site in the 
action (13), so that in the many-body problem the in- 
teracting d— orbital is the only explicit degree of free- 
dom. We then solve the impurity problem using the stan- 
dard Hirsch-Fye QMC algorithm 23 , where the imaginary 
time interval [0, 0\ is discretized in L time-slices of width 
At = f3/L (where /3 is the inverse temperature). We set 
UAt < 1 to limit the systematic errors introduced by 
the Trotter decomposition. The precision of the calcu- 
lations then basically depends on two remaining factors, 
the statistical error and the criterion for the convergence 
of the solution of the DMFT equations. For the former, 
we typically perform 10 sweeps. When required, we may 
do up to 10 6 sweeps, such as near the Mott transition, or 
to compute the analytic continuation of data to the real 
axis using the maximum entropy method 24 . The quality 
of the convergence is controlled by monitoring the be- 
havior with iteration number of Gd d {iuj\), the imaginary 
part of the d— electron Green's function at the first Mat- 
subara frequency, which shows the largest variations. We 
stop the DMFT iterations when the fluctuations of this 
quantity become of the order of the QMC statistical error 
and remains stable for a few more iterations. In generic 
regions of the parameter space we have studied, the solu- 
tion converge in less than 20-30 iterations, but hundreds 
may be necessary close to a phase boundary. 

The ED algorithm is based on the representation of the 
cavity Green's function by finite number auxiliary atomic 
sites 6 . They conform the "bath" or environment of the 
local impurity. In general, an infinite number of sites 
may be required to faithfully represent the dynamic en- 
vironment, however, this is not possible to do in practice. 
Therefore, one has to adopt a strategy to best represent 
the environment with a finite number of auxiliary sites. 
One may use two different "geometries" to represent the 
bath: either the "chain" geometry, as described in Ref. 25 
or the "star" geometry, as described in Refs. 26,27. In 
both cases, the effective impurity problem consists of a 
central impurity site, composed of an explicit d— orbital 
and a p— orbital, where, by virtue of (12) only the latter is 
connected to the bath. In practice, less than 10 sites can 
be dealt with in this method and the ED is performed at 
T — using the Lanczos technique, which is convenient 
to obtain the Green functions. In the case of the chain, 
parameters of the auxiliary atomic sites can be obtained 
in terms of a continued fraction expansion of the com- 
puted Green's functions 6 ' 25 . On the other hand, in the 
case of the star geometry, the parameters are obtained 
by a x 2 minimization of the difference between the com- 
puted local Green function and a finite size parametriza- 
tion of the cavity 6 ' 26 . 

The ED method can be substantially improved by 



6 



supplementing it with the Density Matrix Renormal- 
ization group (DMRG) technique 28 ' 29 ' 30 . Several ED- 
DMRG procedures for the solution of the DMFT equa- 
tions have been proposed recently 28,30 ' 31,32 . The ED- 
DMRG method that we use here is in essence identical 
to the ED with the linear chain 31,33 . Since the linear 
geometry is perfectly adapted for the DMRG procedure, 
we can "grow" the bath to contain a higher number of 
auxiliary sites with respect to the standard ED 30 . In 
practice, we use up to 40 sites. Nevertheless, we point 
out that despite the large number of sites the spectral 
functions on the real frequency axis do not become con- 
tinuous but still show a discrete multipole structure. We 
also observe that the number of poles is roughly similar 
to the number of sites, but the spectral weight remains 
rather concentrated in a relatively small number of poles, 
emphasizing the discrete nature of the Green's functions. 
However, this does not prevent that the agreement on the 
Mastubara imaginary frequency axis is usually excellent 
(e.g. see Fig. 3). 



III. RESULTS 
A. The Mott insulating state 
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FIG. 3: Density of states for the p— and d— electrons (dashed 
and solid line) as obtained from ED-DMRG calculations (with 
a 40 sites chain) at Ao = 1, t p d = 0.9, and ritot ~ 3. Top panel 
has U = 0.5 and /i = 0.612. Bottom panel has U = 2 and 
/i = 1.029. The arrow indicates the width of the Mott gap 
Am- Insets: imaginary part of the p— and d— electron Green's 
functions. Data are from QMC at T = 1/128 (open symbols) 
and ED-DMRG (lines). The finite (zero) value at uj n — > 
shows the metallic (insulating) character of the solution at 
U = 0.5 (U = 2). 



In this section we shall present our results for the Mott- 
Hubbard regime. In the Hubbard model, the system un- 
dergoes a Mott metal-insulator transition when the den- 
sity is n = 1 and the interaction strength U becomes of 
the order of the bandwidth 6 . Here, a similar phenomenon 
is expected as U is increased to a value of the order of 
the effective bandwidth at the Fermi energy ~ tf^/A 
and keeping rid at about 1. As described before, the 
Mott state is realized by setting a value of the interac- 
tion larger than the effective bandwidth U > t^ d / A and 
keeping the total occupation n to t = n p + rid, = 3. Simi- 
larly to the Hubbard model case, in the PAM there is an 
on-site Coulomb interaction acting on d— orbitals. This 
interaction punishes the double occupation of d sites, and 
consequently favors the tendency to localization and to 
magnetic moment formation of the d electrons. 

Notice that we have chosen in the definition of our 
Hamiltonian (1) and (11) the so called non-magnetic form 
for the interaction. This is motivated by the fact that 
we are interested in the Mott physics of a paramagnetic 
correlated state at rid ~ T i-e. = rid[ ~ 1/2. There- 
fore the interaction term U [rid] ~ i) { n di ~ h) approxi- 
mately cancels at the Hartree-Fock level. This allows to 
obtain some immediate physical insight. The cancella- 
tion implies that for low values of U, where first order 
perturbation holds, the interacting density of states of 
the model remains essentially identical as in the non- 
interacting case (see Fig. 2, upper panels). Therefore, 
the position of the correlated narrow band remains ap- 
proximately fixed at the Fermi energy. Since the position 
remains unrenormalizcd, at higher values of U one would 
expect that the narrow band splits, forming a lower Hub- 



bard band and an upper Hubbard band, below and above 
the Fermi energy respectively and both carrying half of 
the spectral intensity of the narrow band. At the Fermi 
energy a large charge gap would then open and the sys- 
tem becomes a Mott insulator. 



1. Opening of the Mott gap 

This scenario is in fact borne out in the actual model 
solution that is shown in Fig. 3. The data corresponds to 
a T = calculation using the ED-DMRG method with 
40 sites in the bath. The values of the interaction are 
U = 0.5 (upper panel) and U = 2 (lower panel). This 
latter value of U is sufficient to drive the system to the 
Mott insulating state. In fact one may estimate that the 
critical U should be about twice the effective bandwidth 
tp d /A (we shall present later on the full phase diagram). 
The insets of the figure contain a comparison of the re- 
sults for the Green's functions in Matsubara frequency 
from ED-DMRG at T = and QMC at the low tem- 
perature T = 1/128. The agreement is very satisfactory. 
The ED-DMRG method also provides the propagators in 
real frequency, we thus plot the more intuitive DOS in 
the main panels of the figure. In the weakly correlated 
case (upper panel), the DOS resembles the non interact- 
ing one: at lower energies —3 % lu < —1.5 there is a band 
with dominant p— character, while at the Fermi energy 
there is a narrower band with mainly d— character. In 
the Mott insulating state shown in the lower panel, the 
DOS consists of three features: similarly as before, there 



FIG. 4: Top panel: rid (solid line), n p (dashed line) and ntot 
(dotted line), as a function of the chemical potential /i, for 
[7 = 2. Data are from QMC calculations at Ao = 1, t p d = 0.9 
and T = 1/64. A large Mott gap opens at U = 2 around 
0-7 < fi < 1.2. Bottom panel: charge compressibility, hi — 
dn/dfi versus fi for the same model parameter as in the top 
panel. 



is a p— like band at high (negative) energies. However, 
the main qualitative difference is that now the narrow 
band at Ep is splittcd in a lower Hubbard band and an 
upper Hubbard band respectively below and above the 
Fermi energy. The Mott-Hubbard character of the tran- 
sition in this parameter regime is seen from the fact that 
both lower and upper Hubbard band have dominant d 
character. Moreover, one also observes that the p com- 
ponent is not negligible, especially in the lower Hubbard 
band. 

We should also mention that the apparent multiple 
peak structure of the main three features appearing in 
the DOS are merely due to the discreteness of the finite 
number of sites used to describe the environment in the 
ED technique. Though we are using as many as 40 aux- 
iliary sites in the environment, the discretization effect 
still remains rather noticeable. Nevertheless, the split- 
ting of the narrow band at the Fermi energy with the 
consequent opening of a large Mott gap Am that signals 
the Mott insulator state is clearly observed. 

The transition from metallic to insulating state with 
the consequent opening of the Mott gap, can be also 
observed from the behavior of the partial p, d and to- 
tal particle occupation n p , rid and n to t respectively, as 
a function of the chemical potential. In Fig. 4 we show 
these quantities and their first derivative k = dn/dfi pro- 
portional to the charge compressibility, for U = 2. The 
plateaux observed in the occupations, with the respec- 
tive vanishing of the compressibility for n to t = and 
ntot = 4, correspond to the completely empty or com- 
pletely full band insulators. The case n to t = 2 corre- 
sponds to the hybridization band insulator, similar to the 
one already discussed in the non-interacting case. How- 



FIG. 5: Top panel: size of the Mott gap Am as a func- 
tion of U, for different position of the p electron band e p = 
—6,-3,-2,-1 (top to bottom). Lower panel: mix-valence 
character v = rid — 1 as a function of U, for the same position 
of the p band as in the upper panel (bottom to top). The 
results are obtained from ED. 

ever, in contrast to the results for the non-interacting 
case, the strong value of the interaction creates additional 
plateaux in the n([i) curves. The new plateaux occur 
when the total number of particles is exactly n to t = 3, 
which is odd and signals the Mott state. 

An important aspect that we should mention is that 
the Mott insulator state occurs where the correlated d 
site has an occupation close to one, but not exactly unity. 
This can be seen is Fig. 4 (top panel), where the Mott 
plateau occurs at rid — 1 + v (and n p — 2 — v), with 
v « 0.13. The specific value of v depends on the hy- 
bridization and therefore one can consider it as a measure 
of the mix- valence character of the Mott insulating state. 
It is the total number of particles exactly equal to three 
(or one hole) what is required for the onset of the Mott 
insulator state. This implies that in this model where 
the mixed p — d valence is explicitly included through the 
hybridization, the Mott localization occurs for a "com- 
posite" object which has a mixed p and d character. 



2. Size of the Mott gap 

The size of the Mott gap is naively expected of be of or- 
der U, since it should mostly reflect the energetic cost to 
doubly occupy the d orbitals. However our results show 
that the size of the Mott gap Am may be substantially 
smaller than the bare value U. In Fig. 5 we plot Am as 
a function of U, for several values of the bare position 
of the p— band e p , which amount to increase the charge 
transfer energy Ao- As the energy of the p— orbitals is 
shifted down to larger (negative) energies, the effective 
bandwidth of the narrow band at the Fermi energy de- 
creases. In addition, the p— electron band becomes es- 
sentially full with rip — > 2 as e p — > — oo (and keeping t p d 
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FIG. 6: Top panel: imaginary part of the d— electron Green's 
functions in the Mott state, from ED-DMRG (for a large finite 
clusters of 30 sites) at T = and fi = 0.729, 0.829, 0.929, 1.029 
(solid, dashed, dotted-dashed, dotted line respectively). Cir- 
cles are the same quantity from QMC at T = 1/64. Lower 
panel: imaginary part of the p— electron Green's functions; all 
model parameters are the same as in the top panel. 



B. Doped Mott insulator 

So far we have shown that the model does have a Mott 
insulator state which similarly as in the Hubbard model 
case develop two incoherent Hubbard bands above and 
below the Fermi energy. However, unlike in the Hubbard 
model case, the size of the gap, i.e. the separation be- 
tween the Hubbard bands, may be substantially smaller 
than U if the hybridization is relatively high. In the 
following section we shall proceed to dope this Mott in- 
sulator with S carriers, with 6 = n to t — 3. As was already 
reported in Ref. 19, we will observe that the insulator to 
metal transitions that can be obtained by either particle, 
5 > 0, or hole doping, 5 < 0, are qualitatively different. 
The former will essentially reproduce the known scenario 
for the Mott MIT that is realized in the DMFT solution 
of one band Hubbard model 6, 33,34,35 . This was to be ex- 
pected since we have tuned the parameters of the model 
to the regime where the identification of the low energy 
physics of the PAM and the one band Hubbard model, 
was expected to hold 40 . However we shall see that, rather 
surprisingly, the hole doping insulator to metal transition 
bears out a qualitatively different scenario. 



fixed). This implies a decrease in the mix- valence char- 
acter of the electrons at the Fermi energy v, as shown 
in the lower panel of Fig. 5. In this limit the size of the 
Mott gap approaches the "bare" value Am ~ U. How- 
ever, it is interesting to observe that the smaller values of 
e p lead to a substantial renormalization of the size of the 
expected Mott gap. This effect can be thought as due to 
an effective screening that the p— electrons provide, or, 
in more naive terms, because the electrons only "feel" 
the repulsive term U during the time they spend on the 
d— orbital, but not when they visit the p site. So as the 
mixed p — d character is increased, the effect of the U is 
renormalized downwards. 

For completeness we show the behavior of the imagi- 
nary part of the Green's functions at low Matsubara fre- 
quency. In Fig. 6 we present numerical results for several 
insulating states obtained varying the chemical potential 
within the Mott plateau. The data were obtained with 
both QMC and ED-DMRG, so also serve to illustrate the 
good agreement between the two methods. Note that the 
imaginary part of both the p and the d electron Green's 
functions goes to zero at u> n — > 0. By analytic continu- 
ation this implies that the p and d electron DOS vanish 
at the Fermi level, consistent with the insulating char- 
acter of the solutions. Accordingly, when the chemical 
potential [i is varied within the Mott plateau, the p— and 
d— electron DOS simply experience a rigid shift in energy. 
However, as we shall see later, there are dramatic changes 
in the DOS lincshapes as the system turns metallic upon 
doping. 



1. particle doping (5 > 0) 

In this section we shall first describe the MIT driven 
by particle doping and demonstrate that it realizes the 
same first order transition scenario as the one in the sin- 
gle band Hubbard model. The metallization of the Mott 
insulator is most directly seen from the changes that take 
place in the density of states. In Fig. 7 we show the evo- 
lution of that quantity for the p— and d— electron com- 
ponents as a function of doping. The data are obtained 
from analytical continuation of QMC results at U = 2 
and low temperature. To compute the analytical con- 
tinuation using the MEM, we used over 10 sweeps, so 
to minimize the uncertainty due to statistical errors. In 
the left top panel we see the DOS for the case where the 
chemical potential [i is close to the upper energy edge of 
the Mott gap, therefore still in the insulator state with 
<5 = 0. The left bottom panel shows the metallic state 
that is obtained when the chemical potential enters the 
upper Hubbard band, doping the parent Mott insulator 
with particles. We observe a broad peak at the Fermi 
energy and a strong transfer of spectral weight from the 
lower to the upper Hubbard band. The d orbital charac- 
ter remains dominant in the DOS upon doping. In the 
left middle panel we show the DOS in the region of small 
doping, very near the transition. The data reveal that 
both the p and d electron components of the DOS show 
a narrow quasiparticle peak at the Fermi energy, flanked 
by the incoherent upper Hubbard band at higher ener- 
gies. This coherent peak carries a small fraction of the 
spectral intensity which is of order i5. In addition, from 
the enhancement of the slope of the self-energy, one ob- 
serves that the quasiparticles acquire a heavy mass. All 
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FIG. 7: Density of states for the p and d electrons (dashed and solid line) at U = 2, Ao = 1, t p d = 0.9. Left: DOS from 
analytically continued QMC data at T = 1/64. Top panel has /i = 1.079, corresponding to a Mott insulating state (5 — 0). 
Middle panel has fi = 1.244, corresponding to small doping 5 — 0.023. Bottom panel has fi = 1.317, corresponding to high 
doping 5 ~ 0.26. Right panels: comparison of DOS from ED-DMRG and QMC data (thick and thin line respectively) at 
6 = 0.264. 



these features are consistent with the MIT scenario found 
in the one band Hubbard model. 

The right panels of Fig. 7 show a comparison of QMC 
and ED-DMRG results for the DOS at high doping. The 
discrete peaks in the ED-DMRG data are due to the fi- 
nite number of sites in the effective bath. Nevertheless, 
the agreement in the distribution of the spectral weight 
of the two methods is very satisfactory. The compari- 
son of the data also allows for a non-trivial benchmark 
of the numerical results and is a useful illustration of the 
advantages and disadvantages of the different numerical 
techniques. ED-DMRG is exact but has the drawback 
of discrete poles structures that persist even for many 
effective atoms in the bath. QMC, in contrast, produce 
smooth spectra but their numerical precision cannot be 
guaranteed due to the uncertainties in the process of an- 
alytic continuation of data to the real frequency axis. 

In Fig. 8 we present the d— electron density of states 
Pd{&) for several values of particle doping. The data show 
the presence of the quasiparticle peak evolving in between 
the two Hubbard bands. One may see that though the 
data show a reasonable systematic evolution some details 
are beyond the precision of the technique, such as the 
precise size of the gap and the form of the bandedges. 

The most substantial confirmation that the particle 
doping driven MIT scenario in the PAM and in the HM 
are in fact qualitatively analogous comes, however, from 
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FIG. 8: Density of states for the d electrons from QMC data 
at U = 2, Ao = 1, t pd = 0.9, T = 1/64. The particle doping 
is 5 = 0, 0.01, 0.02, 0.04, 0.11, 0.21, 0.32 (bottom to top). 



the observation of the hysteresis effect in the particle 
number curve. The hysteresis is a hallmark of the 
first order nature of this doping driven transition and 
it was observed and studied in detail in the Hubbard 
model 34,35 ' 36,37 ' 39 . We also find it here in the PAM, and 
it is most clearly appreciated in the behavior of versus 
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FIG. 9: Hysteresis loop for the rid(fJ>) curve at the upper 
energy edge of the Mott insulating state. Data are obtained 
from QMC for T = 1/64 at U = 2, A = 1, t pd = 0.9. The 
arrows indicate the solution obtained by following the insulat- 
ing (circles) and the metallic solution (squares). Inset: root 
mean square (RMS) of the imaginary part of the d— electron 
Green's function value at the first Matsubara frequency as a 
function of the chemical potential. It shows the enhancement 
of the fluctuation near the two critical values of chemical po- 
tential, where the continued solution ceases to exist. 



\i. There is a strong dependence of the rt^/i) curves as 
the temperature is lowered, signaling strong correlation 
effects been active with very low energies. In Fig. 9 we 
show the occupation of the d electrons as a function of 
the chemical potential ^, obtained from QMC. The main 
panel shows a detail of the QMC data at the low temper- 
ature T = 1/64, where the hysteresis loop can be clearly 
seen. These results were verified using the T = ED 
technique. 

In order to observe the hysteresis cycle, we use as a 
seed for the iterative procedure the converged solution 
from the previous set of parameters 34,36 . Thus, the solu- 
tions can be continuously "followed" in parameter space, 
until it shows a sudden jump. The discontinuous jump 
occur at the approximate position of the spinodal lines 38 . 
Similarly as in previous studies on the Hubbard model 36 , 
the hysteresis defines a region of parameters where two 
solutions of the DMFT equations can be stabilized. The 
true physical transition should occur where the free en- 
ergy of the solutions cross. The precise determination 
of that line is beyond the scope of our present study and 
probably would require further refinement of the numeri- 
cal techniques due to the very low energy scales involved. 

One key point that provides further support to our re- 
sults is the critical slowing down phenomenon observed in 
the QMC calculation at the phase boundaries of the co- 
existence region 34,39 . This phenomenon is characterized 
by an enhancement of the number of iterations required 
to achieve self-consistency and also by an enhancement 
of the statistical Monte Carlo fluctuations that reveal the 
shallowness of the energy landscape when two solutions 




FIG. 10: Main panel: imaginary part of the d— electrons 
Green's function as a function of the Matsubara frequency 
for A = 1, t pd = 0.9, U = 2, T = 1/64 and p = 1.234, 
corresponding to tiny particle doping S — ntot — 3 ~ 0.01. 
Open symbols correspond to the metallic solution, while full 
symbols to the insulating one. Inset: the correspondant 
d— electron density of states pd(u). Thick line is the metallic 
solution, thin line is the insulating one. 



merge. The root mean square deviation of the lowest fre- 
quency component of the d— electron Green's function is 
plotted in the inset of Fig. 9 as a function of \i for the 
low temperature T = 1/64. As we approach the phase 
boundaries of the coexistence region, we can see that the 
root mean square deviation increases. Starting with an 
insulating solution, if we increase the chemical potential, 
the RMS increases until we reach a critical value of /i, 
where the insulating solution disappears (open circles). 
Similarly, starting form a metallic solution at low dop- 
ing, if we reduce the chemical potential, the RMS grows 
until a critical value of /i, where the metallic solution 
disappears (open squares). 

Evidently, the coexistence of solutions can also be ob- 
served from the behavior of other quantities, such as 
the double occupancy, or the low frequency part of the 
Green's functions. The latter is shown in Fig. 10 for both 
p— and d— electron components. 

These results support the claim that the MIT scenario 
for 5 > in the PAM is completely analogous to the one 
found in the Hubbard model investigations 33,34 ' 35 . 



2. hole doping (5 < 0) 

In the previous section we showed that upon particle 
doping, S > 0, the Mott MIT in the PAM realizes the 
same physical scenario as the one observed in the single 
band Hubbard model. While in the latter the particle- 
hole symmetry implies an identical transition for 6 < 0, 
we shall see that this is not the case in the PAM. A key 
point to appreciate is that we shall keep all model param- 
eters fixed, with the obvious exception of the chemical 
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FIG. 11: Density of states for the p— electrons (dashed line) and d— electrons (solid line) for (7 = 2, Ao = 1, t p d = 0.9. Left: 
DOS from analytically continued QMC data at T = 1/64. Top panel has fi — 0.829, corresponding to a Mott insulating state 
(S = 0). Middle panel has fi = 0.529, corresponding to small doping S = ntot — 3 = —0.04. Bottom panel has /i = 0.329, 
corresponding to high doping S — —0.28. The QMC data show the appearance of a broad quasi-particle peak (ZRP) at the 
Fermi level, in addition to the lower and upper Hubbard bands (LHB, UHB) at high energies. Right panels: comparison of 
DOS from ED-DMRG and QMC data (thick and thin line respectively) at the heavy doping 8 = -0.28. The ED-DMRG data 
are obtained using an environment of 30 sites. 



potential, which controls the occupation. Therefore, if 
the chosen model parameters led, for S > 0, to the iden- 
tification of the PAM with the HM physics, then one may 
also expect that this would be the case for S < as well. 
Rather surprisingly this turns not to be the case 19 . In 
this section we shall describe the main physical behavior 
of the model for the hole doping driven MIT, and in the 
next we shall argue about the origin of this unexpected 
result. 

We begin by showing the effect of hole doping in the 
DOS. In Fig. 11 we plot the change in the p and d compo- 
nents of the DOS as the system evolves from Mott insula- 
tor to a hole-doped metallic state (left panels, top to bot- 
tom) . In the insulator state we observe that the chemical 
potential is located within the correlation gap, and the 
lower and upper Hubbard bands can be well appreciated. 
As we already discussed before, in this Mott insulator 
state the DOS at low frequencies has mostly d electron 
character, since the d orbital was initially located at the 
Fermi energy. Upon hole doping, the chemical potential 
moves within the lower Hubbard band. The metalliza- 
tion produces a wide and strong quasiparticle peak at the 
Fermi energy. At small doping, the low frequency part of 
the spectrum has a characteristic three peaks structure: 
the lower Hubbard band around lo ~ —0.4, the quasi- 



particle peak that crosses the Fermi level, and the upper 
Hubbard band at lo ~ 0.8 An interesting aspect to appre- 
ciate is that the quasiparticle peak has, as before, a larger 
d-character, but, in addition, now it also has a substan- 
tial p-electron component. We shall see that this aspect 
will be consistent with the interpretation of the quasi- 
particle peak now emerging not from mere derealization 
of (i-electrons as before, but from the derealization of a 
composite object that involves a p-hole and a d-electron 
spin. 

We can also observe a transfer of spectral weight, with 
an increase of the relative intensity of the lower Hubbard 
band plus quasiparticle, at the expense of a decrease of 
the upper Hubbard band. The structure that is seen 
below lo ?s —1.5 corresponds to the fully filled band which 
remains with predominantly p— electron character. The 
results shown were obtained with high quality QMC data 
(at least 10 6 sweeps to reduce the statistical errors) and 
performed the analytic continuation to the real axis using 
the maximum entropy method 24 . 

In the right panels of Fig. 11 we present a comparison 
of QMC and ED-DMRG results for the DOS at heavy 
doping S = —0.28. Similarly as before (cf. Fig. 7), the 
apparent multi-peak substructures of the ED-DMRG re- 
sults are not physical and only due to the discrete num- 



12 




FIG. 12: Density of states for d electrons for several values of 
hole doping, from analytically continued QMC data at U = 2, 
A = 1, t p d = 0.9, T = 1/64. From bottom to top, the hole 
doping is 8 = 0, -0.01, -0.03, -0.04, -0.09, -0.15, -0.28. 



ber of poles that result from a finite number of atomic 
sites in the auxiliary bath. Compared to the standard 
ED, this method produces much smoother spectra with 
greater detail due to the dramatic increase in the number 
of poles. However, as is also evident from the data, the 
discreteness due the finite size representation of the bath 
remains a shortcoming of the method. The overall com- 
parison of the two methods remains, nevertheless, very 
satisfactory and serves to illustrate their relative advan- 
tages and disadvantages. 

To complete this study, we show in Fig. 12 the de- 
tailed evolution of the quasi-particle peak at the Fermi 
level plus the Hubbard bands as a function of doping. 
The data correspond to the d— electron density of states 
obtained from analytic continuation of QMC data. 

The emergence of a quasiparticle peak at the Fermi 
energy for the hole doped case may seem, at first sight, 
similar to the metallic state obtained from particle dop- 
ing. However, this will not turn to be the case 19 ' 41 and 
the physical origin of the two quasiparticle excitations 
will be shown to be qualitatively different. The reason 
for this unexpected asymmetry will be discussed in the 
next section. 

In order to fully underpin the nature of the order of 
the hole doped transition, we now look for hysteresis ef- 
fects. Thus, as we did before, we continuously follow 
the solutions in parameter space. First, we start from 
the insulator and lower the chemical potential till we ob- 
tain a significantly doped metal; and then, we start from 
the metal and increase /i until we reach again the insu- 
lator. Neither our QMC numerical simulations nor the 
ED-DMRG studies showed any indication of hysteresis 
effects. The QMC data down to T = 1/128 are shown 
in Fig. 13. Compared to the results for particle doping 
(inset), the present ones show a negligible temperature 
dependence. Thus, up to our current numerical capacity 




FIG. 13: Partial occupation of the d electrons, rid, as a func- 
tion of the chemical potential fi at U — 2, Ao = 1, t v d — 0.9. 
Data are obtained from QMC calculations. Main panel has 
T = 1/128. No trace of coexistent solutions are found. In- 
set has T = 1/16, 1/32, T = 1/64 (dotted, dashed and full 
lines). Notice the strong temperature dependence on the par- 
ticle doped side of the transition (where only one branch of 
the hysteresis cycle at T = 1/64 is shown), and the much 
less temperature dependence on the hole doped side of the 
transition. 



we have to conclude that the metal-insulator transition 
in the hole doped case is of second order. As no ev- 
idence of coexistence has been found, neither in finite 
or zero temperature calculations, we are led to the con- 
clusion that the hole doping driven transition at T = 
is continuous, i.e. has a second order character. Thus, 
it is qualitatively different from the particle doped case, 
and, consequently, also qualitatively different from the 
Hubbard model scenario. Of course, we cannot rule out 
the eventual existence of tiny energy scales (which may 
modify our proposed scenario) which remain beyond the 
numerical precision of our methodology. This issue can- 
not be resolved nor by ED or ED-DMRG data due to the 
finite frequency cut-off set by the finite size of the clus- 
ters to diagonalize. Resolving this issue would probably 
require NRG study. However this is not fully clear since 
NRG method requires a good separation of energy scales, 
which is not the case here. 



IV. DISCUSSION 

Our results for the doping driven MIT, that arise from 
either particle or hole doping of the Mott insulator, show 
a qualitative asymmetry, thus questioning the expected 
mapping of the PAM onto the Hubbard model. In this 
section we shall address in more detail this issue from the 
perspective of the physical nature of the two MIT taking 
place in the PAM. 
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FIG. 14: Phase diagram of the PAM in the U — n plane that 
maps the Mott insulator (MI) and its boundaries with the 
metallic state (M). The boundary between the metal and the 
band insulator (BI) at ntot = 4 is also shown. The boundary 
lines are for T = ED star geometry calculations at Ao = 1 
and t p d = 0.9. The dashed line denotes the region of the 
parameters where there is coexistence between two solutions, 
one metallic-like and the one insulating-like. 



A. Phase diagram 

To begin, we map out the phase diagram in U — fx 
parameter space to explore the respective ranges of the 
first and the second order transitions. Our results are 
summarized in Fig. 14, that was constructed using finite 
T quantum Monte Carlo and T = ED data. The latter 
were obtained using the "star-bath" geometry, which we 
found better suited 6 ' 26 than the "chain-bath" geometry 
to study the possible coexistence between the metallic 
and the insulating regions. 

The phase diagram in the U — fi plane shows a cen- 
tral V-shaped Mott insulator region with the correlated 
metallic phases for particle and hole doping, respectively 
to the right and to the left. A threshold value for the 
strength of the interaction U (tip of the V-shape bound- 
ary) is required to obtain a Mott insulator state. This 
threshold depends on the value of the "bandstructure" 
parameters Ao and t p d- This feature is analogous to the 
existence of critical value of the ratio U/D in the one 
band Hubbard model 6 . 

The central V-shaped Mott insulating region shows a 
remarkable asymmetry comparing the hole and the elec- 
tron sides. In contrast, in the one band Hubbard model, 
due to the particle-hole symmetry, the V-shape onset of 
the Mott insulator is symmetric respect to the tip and 
behaves as U c ~ U c i ± 2/i, where U& is the value of 
the interaction at which the insulator disappears 6 . In 
the PAM, the transition line to the particle doped side 
behaves as well like U c ~ const + 2/i. However, the tran- 
sition line to the hole doped metal is almost vertical. 

We also mapped out the order of the transition along 
the boundary lines. Upon hole doping, both QMC cal- 



culations down to the low temperature T = 1/128, and 
ED calculations at T = show no trace of coexistent so- 
lutions along all the transition line, indicating a second 
order transition. On the other hand, upon particle dop- 
ing, the dashed line in the phase diagram displays the 
region of the parameter space where, for sufficient low 
temperatures, an insulating state coexists with a metal- 
lic state. In fact within ED method at T = we find 
that the DMFT equations have two different solutions 
all along the transition line. Within QMC, and down to 
the low temperature T = 1/128, we found a coexistence 
region only for a sufficiently large strength of the repul- 
sive interaction U (approximately U ~ 2 for our specific 
choice of parameters). However, for smaller values of 
U, the phase boundary remains strongly dependent on 
temperature. This suggests that the temperature below 
which there is a coexistence between metallic and insu- 
lating state decreases rapidly approaching the tip of the 
V-shape. This should be expected, since the two bound- 
ary lines, to holes and particle doped metal, eventually 
merge at the tip. Similarly, we also expect the T = 
width of the coexistent region to become narrower as one 
approaches the tip of the V. However, obtaining detailed 
and precise numerical data there turned out to be beyond 
our methods' capabilities. 

We note that at the Mott insulator state (<5 = 0) the 
occupation of the conduction p— electrons is almost sat- 
urated at rip = 2 — v (with v ks 0.1). On the other hand, 
the occupation of the non dispersive d— electrons which 
carry strong magnetic moments due to the on-site repul- 
sion is close to one, with rid = 1 + v (cf Fig. 5). There- 
fore, there is only a small number ~ v of p— holes avail- 
able to screen a number of order one of d magnetic mo- 
ments. Therefore, a natural issue to consider is whether 
the physics associated to the "exhaustion problem" of 
Nozieres 9 ' 42,43 (see also related works on the Kondo lat- 
tice model 44 ' 45 and on the PAM 46 ' 47,48 ) may play a role 
in the different transitions at S greater or smaller than 0. 

An important point to realize is, however, that the 
exhaustion situation is even more extreme on the par- 
ticle side than on the hole- doped side. In fact, while 
the number of d— electron remain always of order one 
on both sides, the number of available p— holes is sub- 
stantially smaller for 6 > with respect to the 6 < 
case. Nevertheless, the Hubbard-like first-order transi- 
tion scenario takes place only on the particle-doped side. 
In other words, the PAM metal-insulator transition sce- 
nario is analogous to the one in the Hubbard model when 
the PAM is even deeper in the exhaustion limit (6 > 0, 
n p w 2). This implies that while exhaustion should play 
a role, it is not obviously responsible for the failure of the 
mapping of the PAM onto the HM for S < 0. 



B. Nature of doped carriers 

We now address the issue of the physical nature of 
the metallic states in this system. Although the low fre- 
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FIG. 15: k = dntot/dfi as a function of doping 8, for sev- 
eral values of the temperatures T = 1/16, 1/32, 1/64 (circles, 
squares, and triangles respectively). Data are obtained from 
QMC calculations at A = 1, t pd = 0.9, U = 2. On the 
particle doped side, we followed the metallic solution. 



quency part of the DOS of both, particle and hole doped 
metals, have a similar three peaks structure with domi- 
nant d character, they realize physically different states. 
For instance, their charge compressibility has a very con- 
trasting temperature dependence. In Fig. 15 we show 
the derivative of the total occupation with respect to the 
chemical potential, k = dn to t/dfi (an observable propor- 
tional to the compressibility) as a function of doping. 
In the Mott insulator (5 = 0) k is zero, indicating the 
incompressible Mott state. As expected, upon doping, 
k increases, indicating that the system becomes com- 
pressible. For small particle doping, the compressibility 
rapidly increases with the temperature, similarly as in 
the one band Hubbard model case. Upon hole doping, k 
is much less dependent on temperature. 

In the PAM the physics of the metallic states is usually 
discussed in terms of the screening between the p conduc- 
tion electrons (or holes) and the magnetic moment of the 
local d— electrons. In fact, at each lattice site a p— and a 
d— electron may form a local singlet, which is the under- 
lying idea in the argument of Zhang and Rice (ZR) in the 
context of high temperature cuprate superconductors 49 . 
However, if on the particle side of the Mott MIT we have 
argued that there are essentially no available holes, then 
the question is, what is screening the d— moments so to 
produce a normal (but heavy) Fermi liquid metal, analo- 
gous to the one in the doped Mott insulator in the Hub- 
bard model? 

The answer to this question is that, similarly to the one 
band Hubbard model case, the d— electrons are screen- 
ing themselves. In fact, the local ZR singlet formation 
does not take place for 6 > simply because there are 
no holes available for screening. Nevertheless, despite the 
high filling of the p— orbitals, the strong hybridization t p d 
still allows for derealization of the d— electrons, through 
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FIG. 16: Schematic representation of the different type of 
carriers that occur upon either (a) particle (5 > 0), or (b) 
hole (S < 0) doping of the Mott state. 



charge fluctuations across the p— sites, with an effective 
amplitude t c s ~ tp d /A. The key physical point in this 
process is that since the p— orbitals are almost full, they 
have a negligible local magnetic moment, so these charge 
fluctuations take place without significant magnetic p—d 
coupling. Therefore the magnetic phase coherence of the 
d— electrons is preserved and, in consequence, a superex- 
change mechanism between neighboring d sites occurs. 
Thus, from the point of view of the d— electrons, they 
have strong magnetic moments, they delocalize keeping 
their quantum mechanical phase through essentially non- 
magnetic p-sites, and therefore also experience antifer- 
romagnctic correlations with nearest neighbor d— sites. 
These physical ingredients arc evidently also realized for 
the carriers in the single band Hubbard model. Therefore 
we can now rationalize the underlying mechanism for the 
mapping of the PAM onto the one band Hubbard model 
at 6 > 0, cf. Fig. 16(a). 

On the other hand, the situation is very different as the 
chemical potential is lowered to dope holes into the sys- 
tem. There, the number of available holes become more 
significant and they can lock with the robust d— magnetic 
moments to form local ZR singlets. However, when the 
d— electrons of these singlets want to delocalize, i.e. hop 
to the neighboring sites and eventually form a quasiparti- 
cle band, they completely loose the information of their 
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spin phase, thus the antiferromagnetic correlations be- 
tween neighboring d do not build up. In consequence, 
Hubbard model like physics does not take place and the 
nature of both the MIT and the ensuing correlated metal- 
lic state becomes fundamentally different. That is the key 
physical reason why the mapping of the PAM with a one 
band HM is no longer valid, cf. Fig. 16(b). 

In order to substantiate the previous qualitative dis- 
cussion, we show in Fig. 17 the d and p local moment 
formation as a function of doping for a large strength of 
the interaction U and for a value slightly above the tip. 
The local moment formation is defined as 

(( m a) 2 ) = (("at ~ n al) 2 ) = n a ~ 2(n a] n ai ) , (17) 

where a — p,d. Notice that the difference between the 
particle occupation n and the moment directly measures 
the double occupancies of the sites. In the Mott insula- 
tor (8 = 0) the local moment of the d electrons, ((m d ) 2 ), 
is large because the d— sites are predominantly single oc- 
cupied due to the effect of U. On the other hand, the 
p— band is almost fully occupied, and thus the local mo- 
ment of the p— electron, ((nip) 2 ), is significantly smaller. 
As one dopes the Mott insulator with particles or holes, 
the relative distribution of the local moments among both 
p and d sites is strikingly asymmetric. Upon particle 
doping, 8 > 0, ({nip) 2 } becomes even smaller, since the 
occupation of the p band get saturated. On the other 
hand, the d moment decreases more rapidly (and lin- 
early) with the doping, since the charge fluctuations be- 
tween the single occupied and double occupied d states 
increase. Therefore these results support the view that 
for the particle doped side of the transition the p— sites 
are magnetically inert and, consequently, the d-electrons 
screen themselves as they form a heavy mass metal. This 
heavy metal paramagnetic state is physically analogous 
to the one realized in the Hubbard model case. For the 
hole doped metal, 8 < 0, the p— electron local moment, 
((nip) 2 ), increases linearly with the doping because holes 
are mostly added to p— sites. In contrast, the opposite 
behavior occurs for the d local moment, which slightly 
linearly decreases with hole doping. Thus, the increase of 
the magnetic character of p— electrons is consistent with 
our argument for the formation of local singlets in the 
hole doped case. 

To fully underpin our hypothesis we compute the mag- 
netic moment correlation between the d— and p— sites, 
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FIG. 17: Top panel: ((m z d ) 2 ) (circles) and {(m^) 2 ) (squares) 
as a function of doping 8 for Ao = 1, t v d = 0.9, [7 = 2 
(full symbols) and U = 1.2 (open symbols). Bottom panel: 
(TripTYid) as a function of doping 8, for the same values of the 
parameters as in he main panel. The results are obtained 
with ED with star geometry. 



Rice singlets in the present model, that only occurs at 
8 < 0. 



C. Formation of Zhang- Rice like singlets 

The origin of the two different MIT scenarios and the 
resulting correlated metallic states can also be argued 
from an energetic point of view. The doping introduces 
new states inside the Mott gap Am, which is renormal- 
ized by the hybridization t p d- These states are a mixture 
of p and d states. An estimate of the energy gain of the 
Zhang Rice singlet formation is 
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This results from the magnetic energy gain from the hop- 
ping of a p— hole (2 — n p = v) on a d— site singly occupied 
(1 — v). Wc should compare the above binding energy 
with the derealization energy of a d— electron between 
two neighboring d— sites, which is proportional to 50 



(m z pm z d ) = ((n d T - n di )(n p] - n pl )) 
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The results are shown in the bottom panel of Fig. 17. 
They illustrate that in fact on the particle doped side the 
magnetic correlations are negligible, however on the hole 
doped side they increasingly grow as the chemical poten- 
tial moves into the lower Hubbard band. The growth of 
the expectation value is commensurate with the increase 
in hole doping and signals that the doped p— holes bind 
magnetically to the local d— electron magnetic moments. 
This represents the formation of the equivalent to Zhang 



This results from the virtual hopping of a d— electron 
to its nearest neighbor d site (i.e. through two p sites). 
A priori this gain is of order one, since the d site are 
approximatively all singly occupied, and 1 — v electrons 
participate in the superexchange process 51 . In fact, in the 
region of parameter we are investigating, the parameters 
t pc [, A and U are of order one. Therefore, the energy 
gain in the derealization of a d— electron (20) is of order 
(l-vf. 
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FIG. 18: Phase diagram of the PAM in the U — /i plane 
that maps the Mott insulator and its boundaries with the 
metallic state. Data are for t p d — 0.9 and Ao = 1,2,3, i.e. 
e p — —1,-2,-3 (solid, dashed and dotted lines). Data are 
obtained with T — ED calculations. 



Upon particle doping, the energy gain of singlet for- 
mation (19) is of order v[\ — i>), thus for small but finite 
v, it is much smaller than the spin exchange energy gain 
(20). Thus we can understand that for particle doping 
the spin exchange energy gain dominates on the energy 
gain of singlet formation. This imply that the metallic 
state realized upon particle doping, is due to the dclocal- 
ization of mostly d electrons. Therefore the d electrons 
play the same role as the single type of carriers in the 
single band Hubbard model and produce qualitatively 
similar MIT scenario. 

On the other hand, upon hole doping, the energy gain 
of Zhang Rice singlet formation can be substantial. In 
fact, in our region of parameters, A is of the same order 
of Am , and the partial cancellation of the denominator 
in (19) explains that for a small but finite v (which is con- 
trolled by the hybridization), the energy gain of singlet 
formation in the hole doped transition case dominates 
on the spin exchange energy gain. Thus we can under- 
stand that, upon hole doping, the ensuing metallic state 
is due to the dclocalization of these composite objects. 
The nature of this metallic state has been studied in de- 
tail in our recent work 41 where we found that it severely 
deviated from the Fermi liquid paradigm. The physical 
reason is that the dclocalizcd holes undergo a strong mag- 
netic scattering from the local magnetic moments as they 
hop from site to site. Thus no coherent behavior for hole 
propagation is realized, at least down to very low tem- 
perature scales, leading to the observed non- Fermi liquid 
character. 



D. Restoring the mapping of PAM onto Hubbard 
model 



We have shown that the new physics in the PAM with 
respect to the Hubbard model comes from the local cou- 
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FIG. 19: Top panel: ((m z d ) 2 ) (circles) and ((nip) 2 ) (squares) 
as a function of doping S for t p d = 0.9 and U approximately 
two times larger than the threshold value to have a Mott gap 
(tip in the V shape of the phase diagram) and e p = —1,-2,-3 
(black, red and green symbols). Explicitly: e p = — 1 and 
U = 2, e p = -2 and U = 0.8, e p = -3 and V = 0.4. Bottom 
panel: (m p m z d ) as a function of doping S, for the same values 
of the parameters as in the top panel. The results are obtained 
with T = ED star geometry calculations. 



pling between the p and d electrons. Thus we expect that 
disfavoring this binding may restore the validity of the 
mapping on the hole doped side of the MIT. To test this 
hypothesis, we lower the energy position of the p— band, 
e p (i.e. we increase the charge transfer energy Ao). We 
already noticed in Sec. Ill A that upon increasing the 
bare charge transfer energy Ao = £d — e p , the Mott gap 
approaches the bare value U and the mix-valence char- 
acter of the electrons is in fact decreased. 

In Fig. 18 we show the phase diagram in the plane 
U — [i for different values of the position of the p band. 
Upon increasing Ao, the threshold value of U to obtain 
the Mott insulator region becomes smaller. This results 
from the fact that the interaction U is competing with a 
decreasing bandwidth ~ tiz d /A. In addition, as expected, 
the V-shaped boundaries of the Mott insulator become 
more symmetric at larger Ao. 

To verify whether the character of the transition can be 
modified on the hole doped side, we observe the temper- 
ature behavior of the particle number, n. Upon increas- 
ing Ao, we find that n versus \i curves becomes more 
temperature dependent, which is a first indication of the 
possible realization of the Hubbard model scenario also 
in the hole doped side. However, due to the reduction of 



the effective bandwidth t 



eff 



tp d /A, the temperature 



below which we may observe the hysteresis cycle in the n 
versus \x curves should be extremely low. With our cur- 
rent numerical capabilities we were only able to obtain 
evidence of a small hysteresis at Ao = 3. 

To complete this study we also computed the local 
moment of p— and d— electrons upon increasing Ao. In 
Fig. 19 we plot the moments (top panel) and the mag- 
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netic binding between the p— and d— electrons (bottom 
panel) as a function of doping for several values of the 
position of the p— band. In the Mott insulator (<5 = 0), 
upon decreasing the position of the p band, {{m p ) 2 ) is 
decreased and ((m^) 2 ) is increased. This is because 
the mix-valence character of the electrons is reduced, 
since the p band become essentially full with n p — > 2 
as e p — > — oo. Thus the doping has just a small effect 
on {(m z p ) 2 ). On the other hand, upon hole doping, the d 
local moments decreases more rapidly when the p— band 
is deep in energy. As a result, the magnetic correlation 
between the p— and d— electrons, (m p m^), shown in the 
bottom panel of Fig. 19, is expectedly punished by higher 
values of Ao . This is fully consistent with the mentioned 
recovery of the mapping of the PAM onto the Hubbard 
model. Therefore we can understand that for lowering 
in energy the position of the p— band the mapping of the 
PAM onto the HM should hold. 



V. CONCLUSIONS 

In this paper, using DMFT, we considered the doping- 
driven Mott transitions of the periodic Anderson model 
when it is set in the Mott-Hubbard parameter regime. 
We discussed the transitions with respect to reference 
case of the well understood scenario realized in the single 
band Hubbard model. 

In contrast to the latter, the PAM has a qualitatively 
different metal-insulator transition for particle or hole 
doping. Upon particle doping of the Mott insulator, the 
metallic state is reached through a first order transition, 
that is analogous to that of the Hubbard model. How- 
ever, upon hole doping the Mott insulator, there is a 
continuous (i.e. second order) insulator-metal transition 
through a quantum critical line in the parameter space 
U-fi. 

We argued that the hole doped metal has delocalized 
Zhang Rice singlets that fail to build substantial superex- 
changc as compared to the Hubbard model (and particle 



doped case). In fact, we discussed the qualitative differ- 
ences between these two transitions, showing that it is 
not due to the physics of the "exhaustion" , but indeed 
is related to the magnetic interaction that develops be- 
tween the two species of electrons in the model. Our re- 
sults on the magnetic correlation between the d— and the 
p— electrons (see lower panel of Fig. 17) show that upon 
particle doping the p— electrons permit the charge fluctu- 
ations and the derealization of the d— electrons without 
magnetic p — d coupling. On the hole doped case, in 
contrast, the system favors the formation of singlet pairs 
p — d. 

Upon increasing the charge transfer energy, we could 
recover the mapping of the PAM to the Hubbard model 
for the hole doped case. This signifies that a substan- 
tial mix-valence character was the key ingredient for the 
realization of the second order transition in this model. 

Our findings may be important looking at the present 
effort to apply DMFT calculation in regard to real 
materials 52,53 . Those studies usually carry the implicit 
assumption of the Hubbard model as the underlying low 
energy Hamiltonian of complex systems. Our work in- 
dicates that the Hubbard model scenario may be ques- 
tionable when the hybridization of the correlated band 
with another band is high. In particular, our work is rel- 
evant for the analysis of the metal-insulator transitions 
of transition-metal oxides, that usually have oxygen or- 
bital mediating the derealization of the d correlated elec- 
trons of the transition-metal. Therefore the role of the 
oxygens band and their hybridization with the localised 
band should be explicitly considered in the investigation 
of the Mott transition. 
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